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PACS 87 . 23 . Ge - Dynamics of social systems 
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Abstract. - For models whose evolution takes place on a network it is often necessary to augment 
the mean-field approach by considering explicitly the degree dependence of average quantities 
(heterogeneous mean- field). Here we introduce the degree dependence in the pair approximation 
(heterogeneous pair approximation) for analyzing voter models on uncorrelated networks. This 
approach gives an essentially exact description of the dynamics, correcting some inaccurate results 
of previous approaches. The heterogeneous pair approximation introduced here can be applied in 
full generality to many other processes on complex networks. 



Introduction. — The voter model is one of the sim- 
plest possible models of nonequilibrium dynamics, as wit- 
nessed by the number of different contexts where it has 
been considered, ranging from probability to theoretical 
ecology and heterogeneous catalysis [l}{3]- 

Each node of a graph is endowed with a binary variable 
(a spin) s — ±1. At each time step a node and one of its 
neighbors are chosen at random and the first node becomes 
equal to the second. This microscopic dynamics gives rise 
to a nontrivial macroscopic phenomenology [3] with the 
formation of correlated domains that tend to grow in time. 
In finite systems this eventually leads to one of the two 
possible ordered states (consensus) with either all spins 
up (+1) or down (—1). In the thermodynamic limit full 
order is reached only in euclidean dimension d < 2. 

From the point of view of statistical physics voter dy- 
namics has the very interesting feature that its overly sim- 
ple rules allow for analytical treatment in many cases. On 
euclidean lattices many exact results are available [HIS]. 
In the past few years, in the framework of the huge in- 
terest around complex networks and dynamics on top of 
them [SHS] , the effect of nontrivial topologies on voter dy- 
namics has been explored [9Ul8] . One interesting effect of 
the presence of a disordered connectivity pattern is that 
slightly different definitions of the model, perfectly equiv- 
alent on lattices, become inequivalent on networks. In the 
direct voter case, one picks up a vertex at random and 



this imitates a randomly chosen neighbor. In the reverse 
voter dynamics [13j or invasion process the order of selec- 
tion is opposite: it is the neighbor that imitates the node 
selected first. An intermediate case is the link-update dy- 
namics [12| , where one has to choose at random first a link 
and then which of the two nodes at the end of the link imi- 
tates the other. Only the link-update rule preserves a cru- 
cial property of voter dynamics on lattices, conservation 
on average of magnetization [T^ , whereas other quantities 
are conserved in the two other cases [TTJ[T3]. On regu- 
lar geometries the three types of dynamics coincide, but 
when the topology is strongly heterogeneous important 
differences arise: on scale-free networks the scaling of 
average consensus time T with the system size N changes 
depending on which of the three versions is considered. 

Theoretical progress on the understanding of these is- 
sues has been made possible by the application of a hetero- 
geneous mean- field approach [TT 1 IT3 l fT6 ] . The basic idea is 
to define the densities at of up spins restricted to nodes of 
degree k and to write down evolution equations for them. 
In this way it is possible to predict that for the direct voter 
the average time T to reach consensus grows sublinearly 
as a function of the system size N if the exponent 7 of 
the tail of the degree distribution is smaller than 3, while 
it is linear for 7 > 3. Instead, for the reverse and link- 
update versions of the dynamics the growth is linear for 
any 7 > 2. Very recently, the theoretical picture for the 
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voter model has been further refined by the apphcation 
of an improved mean-field approach, a homogeneous pair 
approximation |20| . taking into account dynamical corre- 
lations among nearest neighbors. This approach allows to 
follow in detail the temporal evolution of the system and 
to compute with remarkable accuracy even the prefactor 
of the law that expresses the average consensus time as a 
power of N. 

Despite these achievements, some open problems re- 
main, calling for further improvement. The first has to 
do with the prediction of a transition for average degree /x 
equal to 2. According to Ref. PD|, while for > 2 consen- 
sus is attained via a quasi-stationary intermediate state 
over a long temporal scale, when /i < 2 order is reached 
exponentially fast. Contrary to this prediction, numerical 
simulations on an Erdos-Renyi graph do not display any 
singular behavior in correspondence oi ^ — 2 20 . An- 
other issue that requires a clarification has to do with the 
use of the same approach for reverse voter dynamics. As 
it will be explained below, the application is straightfor- 
ward, but the agreement with numerical simulations is not 
complete. 

In this paper we go beyond the homogeneous pair ap- 
proximation introduced in Ref. |20| by allowing correla- 
tions between nearest neighbor nodes to depend on their 
degree (heterogeneous pair approximation). In practice 
we lift the assumption of Ref. |25 that the probability for 
a link between two nodes to be active (i.e. to connect 
nodes in different spin state) does not depend on the de- 
gree of the nodes. This is analogous to the introduction 
of the heterogeneous mean-field approach [8l[2T| necessary 
for dealing with the SIS model (and many others) on scale- 
free networks. With the heterogeneous pair approximation 
we find, for the direct voter dynamics, results in general 
very close to those obtained with the homogeneous pair 
approximation but we get rid of the spurious transition 
at /i = 2, in agreement with computer simulations. Also 
for the reverse voter and link-update dynamics we obtain 
results in very good agreement with numerics, showing 
that the heterogeneous pair approximation fully captures 
the behavior of voter models. These successful results call 
for the application of the heterogeneous pair approach to 
other nontrivial dynamical processes on networks, includ- 
ing also epidemiological models or evolutionary games [22, . 
While the detailed implementation will differ depending 
on the case, the idea is completely general and may lead 
to substantial progress with respect to the heterogeneous 
mean-field approach. 

Problems with the homogeneous pair approxi- 
mation. — The homogeneous pair approximation intro- 
duced in Ref. [20 is based on the derivation of the equation 
of motion for the quantity p, the fraction of active links, 
i. e. edges connecting nodes in opposite spin state. In the 
derivation, it is assumed that the probability for a link in 
the system to be active does not depend on the degree of 
the nodes it connects, and is equal to p. Under this hy- 



pothesis the equation of motion for p is derived, involving 
the magnetization m and the first moment of the degree 
distribution /i. As mentioned before, the behavior pre- 
dicted changes dramatically as a function of p. For /i < 2, 
p goes to zero exponentially fast in time, indicating that a 
fully ordered configuration (consensus) is quickly reached: 
the average time T to reach consensus does not depend on 
the system size N so that T/N 0. For p > 2 instead it 
is predicted that, over a very short temporal scale, a quasi- 
stationary state is reached, characterized by a density of 
active links 
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The superscript S indicates, throughout the Letter, quan- 
tities in the quasi-stationary state. The eventual behavior 
of the system consists in the erratic variation of m and 
p^ that jointly fluctuate in time while obeying Eq. ([1]), 
until a fluctuation leads to the absorbing state \m\ = 1, 
p — 0, corresponding to full consensus. Other dynamical 
observables, describing the ordering process, are expressed 
using Eq. ([T]), as for example the average consensus time 
that, for uncorrelated initial conditions with magnetiza- 
tion m = 0, is 



T = 



p^N 



2pS{m 0)p; 



■ln(2) 



{p - 2)p2 



ln(2), (2) 



where p2 is the second moment of the degree distribution. 
While Eq. ([T|) agrees with remarkable precision with the 
outcome of numerical simulations for large values of p, 
significant discrepancies arise on many types of networks 
when p gets close to 2. For example, considering as a sub- 
strate an Erdos-Renyi graph, the divergence of T/N for 
p ^ 2 predicted by Eq. ([2]) is not seen in numerical simu- 
lations [20] (Fig. [1]). This indicates that, despite the very 
good performance for large values of p, the homogeneous 
pair approximation of Ref. [20] qualitatively fails in the 
limit p ^ 2. 

The same indication emerges when the approach is 
applied to the reverse voter dynamics. For the quasi- 
stationary value of the density of active links it is straight- 
forward to obtain 
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This formula coincides, as expected, with the result for the 
direct voter model [Eq. ([1])] when the degree distribution 
is a delta- function so that p2 = p'^- However, a care- 
ful comparison with numerics shows that Eq. ([3]) is not 
correct when p2 ^ p^ (see Figure [5]). The parabolic de- 
pendence on magnetization remains (not shown), but the 
value for TO = predicted by Eq. ([3]) (diamond symbols) 
is considerably larger than what simulations give (circles) , 
unless 7 is very large. Here, as in the rest of the paper 
unless specified otherwise, we compare the results of the 
theoretical treatment with simulations performed on the 
uncorrelated configuration model (UCM) [SJ , that gener- 
ates graphs with power-law degree distribution k~'^ 
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Fig. 1: Average time to reach full consensus for direct voter 
dynamics on an Erdos-Renyi graph, as a function of the av- 
erage degree fi, starting from a symmetric uncorrelated ini- 
tial configuration ak{t = 0) = 1/2. The curve for the het- 
erogeneous mean-field approach is obtained using the formula 
T/N = /i^//i2ln(2) from Ref. [11]. The curve for the homo- 
geneous pair approximation is obtained using Eq. ([2}. The 
curve for the heterogeneous pair approximation is obtained us- 
ing Eq. pOI) . For finite and fi very close to 1, the network 
is separated in many nonextensive disconnected components. 
This is the origin of the discrepancy between numerical results 
and the theoretical prediction in an interval that shrinks to 
zero as the system size diverges. Each point is the average 
over 1000 realizations. 
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Fig. 2: Reverse voter dynamics: quasi-stationary value {rn — 
0) as a function of the system size for three different values of 
7 (top to bottom, 7 = 2.5, 7 = 4, 7 = 8). 
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Fig. 3: Plot of the probability that a link connected to a node 
of degree k is active in the quasi-stationary state (and for m = 
0) for reverse voter dynamics taking place on an uncorrelated 
network with 7 = 2.5 and size iV = 10*^. The black symbols are 
the results of a numerical simulation. The dashed red line is 
the prediction of the homogeneous pair approximation, Eq. Q, 
the blue solid line is the prediction of the heterogeneous pair 
approximation. 



and upper bound of the degree distribution equal to iV^/^, 
so that automatically there are no self-loops nor double- 
links . The lower-cutoff of the degree distribution is 3. 

The origin of the discrepancies in Figs. (H]) and ([2]) is not 
hard to identify. It is sufficient to evaluate numerically the 
quantity , defined as the fraction of neighbors of a node 
of degree k that are in a different spin state in the quasi- 
stationary state, and plot it as a function of k (Fig. [3]). 
While the approach of Ref. [5D] assumes this quantity to 
be the same for all nodes, a substantial dependence on k 
shows up. Interestingly, a numerical check indicates that 
also for direct voter dynamics has a sizeable depen- 
dence on the degree (Fig. |4]) . Based on this evidence, we 
now generalize the pair approximation approach by con- 
sidering pk.k', i-e. allowing the density of active links to 
depend explicitly on the degree on the connected nodes. 
From pk^k' it is possible to determine all other quantities 
of interest. The density pk of active links connected to a 
node of degree k is pk = Y^f,, Qk'Pk,k', where Qk = Pkk/p 
is the degree distribution of the neighbours of a randomly 
chosen node. The total density of active links p is likewise 
given by p = J^k QkPk = Y.kM QkQk'PkM- 

Heterogeneous pair approximation for the direct 
voter model. — The first goal is to determine the equa- 
tion of motion for pk,k' ■ This quantity is modified if the 
flipping node has degree k and one of its neighbours has 
degree k' (or vice versa). Let us assume that the flipping 
(first selected) node has degree k and call k" the degree of 
the copied (second selected) node. It is useful to consider 
separately the two cases where k" ^ k' or k" = fc'. 
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Fig. 4: Plot of the probability that a link connected to a node 
of degree k is active in the quasi-stationary state (and for m = 
0) for direct voter dynamics taking place on an uncorrelated 
network with 7 = 2.5 and size A'^ = 10®. The black symbols are 
the results of a numerical simulation. The dashed red line is 
the prediction of the homogeneous pair approximation, Eq. {T]), 
the blue solid line is the prediction of the heterogeneous pair 
approximation. 



In the first case the variation Apk^k' for a single dynami- 
cal step (occurring over a time At = 1 /N) is determined as 
follows. The probability that a node in state s and degree 
k flips is given by the probability Pk that the first node 
selected has degree k times the probability (t{s) that it is 
in state s, times the probability Qk" that the second has 
degree k" multiplied by the probability pk.k" /[2cr{s)] that 
the link connecting the two is active. One has then to mul- 
tiply this quantity by the associated variation of the frac- 
tion of active links between k and k' (Fig. [S|). Among the 
k — 1 other links of the flipping node, the number of those 
connecting to a node of degree k' will be j distributed ac- 
cording to a binomial i?(j, fc — 1) with probability of the 
single event equal to Qk'- In their turn, only n out of these 
j links will be active, with n binomially distributed with 
single event probability pk.k' /[^(^(s)]. Finally one has to 
multiply by the variation of pk^k' when n out of j links go 
from active to inactive as a consequence of the flipping of 
the node in k. This is given by the variation of the number 
of active links [{j — n) — n] divided the total number of 
links between nodes of degree k and fc' {NfxQkQ'i^). One 
has then to sum over fc" ^ k' , s, j and n, obtaining 



Apk,k' = Pfc^CT(s) ^ Qk' 



k"^k' 



Pk,k" 

2a{s) 



(4) 
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By performing explicitly the summations (and using 
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Fig. 5: Direct voter dynamics: Illustrative representation of a 
configuration contributing to dpsj/dt. The first node selected 
has degree k — 7 and the second has degree k" = 4, different 
from k' — 3. Filled (empty) circles denote spins in state s — +1 
(s — —1). The variation to the density pk,k' is due to the 
J = 3 links that connect the flipping node to nodes of degree 
fc'. Before the flipping of the flrst node selected, n = 1 of the 
links are active and j — n = 2 are inactive. After the flipping 
the state of links is reversed and the variation in the number 
of active links between nodes of degree fc and fc' is j — 2n = 1. 

J2s = 4/(1 — TO^)) the formula becomes 



Apk,k 
At 



-= E Q 

k"^k' 



k"Pk,k" 



k 



1 



TPk,k' 



(5) 

When fc" = fc', the value of Apk,k' is similar to Eq. (j4]) 
with (obviously) Qk> instead of Qk", no sum over fc", and 
in the numerator of the last factor j + l — {n + l) — {n + l) = 
j — 2n — l, because there are j -1-1 links to nodes of degree fc, 
n -|- 1 of which are active in the initial state and inactive in 
the final. Summing up the two contributions and adding 
the symmetric terms with fc and fc' swapped, we get 
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dt 
2pk 
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= Pk- 
k - 1 
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From this equation, by summing over fc and fc' one obtains 
the equation of motion for the total density p 



Tt T.Q^Q^ 



k,k' 



dpk,k' 

dt 
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P 
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P2 



(7) 



where tt ^ J2k PkPk, tt* = ^^fc PkPl and pi = Y.k QkPl- 
It is easy to see that, if one assumes pk = p, i.e. that the 
density does not depend on the degree (and as a conse- 
quence TT = p, P2 = ■"'2 = P^); Eq. ([7]) coincides with the 
analogous equation of Ref . [20] . 

After a time scale of order unity the quasi-stationary 
state is established. From dpk,k'/dt = we obtain 



Pk,k' 



Pk k 



nS k^ 
Pk' k' 



s fc-i 
fe 



■pi'^ 



(8) 



p-4 



Heterogeneous pair approximation for voter models on networks 



From this equation it is evident that ah pf are pro- 
portional to 1 — m^. Imposing the consistency condition 
Pk = J2k' Qk'Pk,k', we have a set of coupled equations for 
for all k. The numerical iterative solution of this set of 
equations allows the determination of the quasi-stationary 
values of all and hence of all pf f., . 

Fig. |4] shows that the heterogeneous pair approximation 
captures extremely well the detailed dependence of pf on 
k. Also the comparison of the quasi-stationary value of 
(recovered by summing Eq. ^ over k and k') with simu- 
lations for several values of 7 and N indicates an excellent 
agreement (not shown). 

The knowledge of the pf allows to determine also the 
consensus time T to reach a fully ordered configuration. 
Following the formalism of Ref. , using the backward 
Kolmogorov equation, one simply needs to express the 
transition probabilities Rk and Lk in terms of the pk- 
In the quasi-stationary state the probability i?f that the 
number of nodes of degree k in state s = increases by 
1 is 



is slightly modified 



= Pk{i - Ok) J2 Qk"Pik"im - <yk)] = p, 



Pk 



(9) 



Plk' 



^l{k-l)+^l,{k' -I) 



[vrf (fc - 1) + 7rf,(fc' - 1)] + p./k + p/k' 



(12) 



where tt/j = Pk' Pk.k' ■ Much in the same way as for 
the direct voter dynamics a set of consistency equations 
for TT^ is solved iteratively. The values of tt^, inserted 
into Eq. (1121) provide p^ ^, , from which in turn p^ and p^ 
can be computed. Fig. |3] shows that the agreement with 
the values of the pf obtained numerically is very good. 
Moreover, Fig. [2] shows that the quasi-stationary value 
of p^ is in fully satisfactory agreement with numerics for 
several values of 7 and N . 

Also in this case the consensus time is computed us- 
ing the backward Kolmogorov equation for T |16j . The 
transition probabilities in the quasi-stationary state are 
obtained from those of the direct case with the replace- 
ment PkQk" Pk"Qk- 



Ru — 



Qk^ 



(13) 



Li. — Ri. . 



The probability Lk that the number is reduced is the same, 
= Rk- Using these expressions one obtains 



T = -t{N) [(1 - u) ln(l - w) - w Inc 



where uj 
and 



J2k Qk<yk{t = 0) is conserved by the dynamics 



t{N) 



2EkPkk'plim^0) 



(11) 



sets the temporal scale over which consensus is reached. 
Notice that with the mean-field assumption, pf (m = 0) = 
1/2, Eq. pi^ coincides with the result of Sood and Redner 
t{N) = N p? / p2- With the homogeneous pair approxima- 
tion pf (to = 0) = (p — 2)/[2(/Lt — 1)] it returns the result of 
Vazquez and Egufluz With Eq. pUj) . valid for any de- 
gree distribution, we have computed the value of the ratio 
T /N for an Erdos-Renyi graph of average degree p. The 
plot in Fig.[T]presents the comparison of this value with the 
prediction of previous approaches and with the outcome of 
numerical simulations. The singular behavior for p = 2 is 
an artifact of the homogeneous pair approximation, that 
is removed by the heterogeneous pair approximation, lead- 
ing instead to a smooth behavior around p — 2 and a very 
good agreement between theory and simulations. 

Heterogeneous pair approximation for the re- 
verse and link-update voter dynamics. The ap- 
plication of the heterogeneous pair approximation to the 
reverse voter model proceeds along the same lines of the 
direct dynamics. The only difference is that the factor 
PkQk" is replaced by Pk"Qk, because the role of the 
flipping and the copied nodes is swapped. As a conse- 
quence, the formula for the quasi-stationary value of pk^k' 



Using these quantities and the fact [131 [16] that the dy- 
namics conserves the quantity uj-i = '^f, Pkcrkk~^ / p-i, 
where p-i —J^k^k/k, one obtains 



(10) T = -r(A^) [tj_iln(a;_i) + (1 - a;_i)ln(l -w-i)] , (14) 



and the temporal scale is 
t{N) 



Nppli 



(15) 



1/2 



For uncorrelated initial conditions ak{t — 0) 
w-i = 1/2, 

T = T(7V)ln(2). (16) 

The comparison of the predictions of Eq. (jl6p with numer- 
ical results is displayed in Fig. [S] and again confirms the 
correctness of the theoretical treatment. 

For link-update dynamics, all calculations proceed ex- 
actly as in the direct voter case, now with the replacement 
of PkQk" with QkQk"- It is straightforward to obtain 



Pk,k' 



pUk-l)+pi'ik'-l) 
j^[pUk-l) + plik'-l)]+2- 



(17) 



The consistency condition pk — X]fe' Qk' Pk.k' is solved it- 
eratively. Once again, the values of pf and pf j,/ obtained 
in this way are in very good agreement with numerics (not 
shown). The transition probabilities for the Kolmogorov 
backward equation are R^ ^ ^ QkP^/2, that yield, 
for the consensus time 



T = t{N) [wo In(wo) + (1 - wo) ln(l - ujq)] , 



(18) 



where the conserved quantity ujq is a{s — +1), the total 
density of nodes in state s — +1 and t{N) — N/[2p^{m = 
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Fig. 6: Reverse voter dynamics: Plot of the consensus time T 
as a function of the system size TV for two values of gamma. 
Filled symbols are for numerical results, solid lines are the the- 
oretical predictions of Eq. p6[) . dashed lines are the mean-field 
predictions of Ref . [TJj . 



0)]. Notice that using p^[m — Q) — 1/2 one recovers 
the mean-field result t{N) — N, i.e. the independence of 
consensus time of network features |13) . However, since 
[m = 0) actually depends on the network, this inde- 
pendence is just an artifact of the mean-field approach. 

Conclusions. — In this Letter, we have introduced a 
heterogeneous pair approximation to deal with the order- 
ing dynamics of voter models on networks. This approach 
allows to determine, with great precision, not only the 
scaling but also the prefactors of the relevant dynamical 
observables (quasi-stationary density of active links, con- 
sensus time) . The heterogeneous pair approximation thus 
constitutes a very good approximate theory for voter mod- 
els on uncorrelated networks. The only point that remains 
to be understood is the origin of the very small discrep- 
ancy (at most of the order of 1%) that still exists between 
numerical and analytical results. 

While the treatment presented here is valid for un- 
correlated networks, generalization to correlated nets is 
straightforward. In particular, if P{k'\k) is the probabil- 
ity that a neighbor of a randomly chosen node of degree 
k has degree fc', the stationarity conditions for the three 
types of dynamics, Eqs. ([8|), (fT2|) and (fT7| . remain the 
same, provided P{k'\k) replaces its expression Qk vahd 
in the uncorrelated case, in the definitions of pk and tt^: 
Pk = Ek' P{k'\k)pk,k-, TTfc = Efc' t^P{k'\k)/k'pkM'- While 
the generalization of the analytical approach is immediate, 
an interesting question deserving further investigations is 
whether the agreement between theory and numerics re- 
mains good even in the correlated case. 

The idea behind the heterogeneous pair approximation 
introduced here is fully general and its application to other 
dynamical models evolving on complex networks is indeed 



a promising avenue for future work. In particular, its role 
could be crucial for understanding cases where heteroge- 
neous mean-field approaches do not work |24] . 



We thank A. Barrat and R. Pastor-Satorras for a critical 
reading of the manuscript. 
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